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The electronic tunneling properties of the two stable forms of trilayer graphene (TLG), rhombo- 
hedral ABC and Bernal ABA, are examined for pn and pnp junctions as realized by using a single 
gate (SG) or a double gate (DG). For the rhombohedral form, due to the chirality of the electrons, 
the Klein paradox is found at normal incidence for SG devices while at high energy interband scat- 
tering between additional propagation modes can occur. The electrons in Bernal ABA TLG can 
have a monolayer- or bilayer-like character when incident on a SG device. Using a DG however 
both propagation modes will couple by breaking the mirror symmetry of the system which induces 
intermode scattering and resonances that depend on the width of the DG pnp junction. For ABC 
TLG the DG opens up a band gap which suppresses Klein tunneling. The DG induces also an 
unexpected asymmetry in the tunneling angle for single valley electrons. 

PACS numbers: 72.80.Vp, 73. 21. Ac, 73.23.Ad 



I. INTRODUCTION 

The discovery of a one atom thick layer of carbon 
atoms, graphene, opened up an entire new field in the 
condensed matter worlds. The electronic properties of 
graphene have an intriguing analogy with ultrarelativis- 
tic particles such as the linear electronic spectrum and 
the occurrence of Klein tunneling 2 - which both have been 
experimentally verified- - -. This analogy leads to the in- 
troduction of the concept of pseudospin 2 giving the carri- 
ers a chirality which is closely linked to these phenomena. 

Graphene multilayers though being bound by a weak 
Van der Waals force posses energy spectra that are fun- 
damentally different from the monolayer case^ 7 -. The low 
energy spectrum is no longer linear and different modes 
of propagation become possible due to the presence of 
multiple energy bands that appear as a consequence of 
the increasing number of atoms in the unit cell. It was 
shown that the behavior of these bands&S, their response 
to an external applied gate voltage^ as well as the trans- 
port properties of the syste m 11 ' 12 strongly depend on the 
way the graphene sheets are stacked. 

The transport properties of bilayer graphene (BLG), 
the thinnest multilayer structure, has been extensively 
studied 2 - ' 13 ' 14 . In contrast to monolayer graphene 
(MLG), it was shown that, within the two band approxi- 
mation, Klein tunneling does not occur, i.e. the transmis- 
sion for normal incidence is practically zero because, due 
to pseudospinorial arguments, the negative energy states 
are cloaked from the positive ones^. However, when 
the full-band Hamiltonian model is considered, it was re- 
cently reported^ that normal transmission becomes pos- 
sible for high enough potential barriers, which makes the 
higher energy bands available for conduction. 

Recently, there has been a growing interest in the study 
of the electronic properties of trilayer graphene (TLG), 



since it constitutes the simplest multilayer system where 
both types of stacking order, Bernal (ABA) and rhom- 
bohedral (ABC) are possible. These stacking types are 
schematically shown in Figs. Ufa) and (b), respectively. 
The mirror symmetric Bernal stacking is the most com- 
mon and can be exfoliated from natural graphite since 
it shares its crystalline structure^ 7 -. Although the ABC 
stacking is less common, it has recently been reported 
that 16% of the synthesized graphite^ and around 15% 
of exfoliated TL G 11 ' 18 has rhombohedral stacking. The 
way in which the layers are stacked influences strongly 
the energy spectrum of the system. It was demonstrated 
that for trilayer systems the occurrence of Klein tunnel- 
ing depends on the staking order, being present only in 
ABC stacked trilayers 1 ^— . Recent experiments have in- 
vestigated the electronic transport in trilayer o 22 ' 23 . 

The tunneling problem of charge carriers in TLG is 
technically more complex which is the reason why only 
very recently this problem was tackled. Kumar et aWs 
calculated the low energy tunneling through pnp junc- 
tions. Unfortunately, the numerical results were shown 
not to be correct^. Here we extend this work to the more 
difficult regime of high energy and high potential when 
several propagation modes are present and thus multiple 
tunneling and reflection channels have to be taken into 
account. We consider both ABC and ABA TLG that are 
affected by a single gate (SG) that is able to locally vary 
the potential on each layer with the same value. Ad- 
ditionally, we calculate the multiband tunneling in the 
presence of a nanostructured double gate (DG) that af- 
fects the potential on each layer separately. We will show 
that this introduces interband scattering even at low en- 
ergy in ABA TLG. 

The paper is organized as follows. In Sec. |H] we in- 
troduce the full band continuum model used to describe 
ABC and ABA TLG and define the different potential 
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FIG. 1: (Colour online) (Top) Two different crystallographic 
structures showing the relative position of the sublattices a; 
and /3i for (a) ABC TLG and (b) ABA TLG. The interlayer 
hopping that is considered in this study is indicated by the 
yellow lines between the yellow marked atoms. (Bottom) En- 
ergy spectrum of (c) ABC TLG and (d) ABA TLG shown 
by the blue solid curves. The red dashed curves in (c) cor- 
respond to the spectrum of the two band Hamiltonian. The 
red dashed curves in (d) correspond to the spectrum of ABA 
TLG with interlayer bias 5 — O.57I. 



profiles that are used in Sec. IIIII In Sec. [TV] we present 
in detail the formalism used to calculate the transmis- 
sion and reflection probabilities for these multichannel 
systems. In Sec. |V]we show the results for transmission, 
reflection and conductance for both stacking sequences 
and considered different potential profiles. Finally, in 
Sec. IVII we summarize our main conclusions. 



which is given by 
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where^ 71 = 377meV is the interlayer hopping param- 
eter. This Hamiltonian is written in the basis of the 
atomic orbital eigenfunctions 



(3) 



where the indices indicate the sublattice associated with 
the respective eigenf unction. The energy spectrum of 
this Hamiltonian is depicted in Fig. [lfc) by solid curves. 
The spectrum consists of six energy bands of which two 
touch each other at k = 0. The two touching bands can 
be described approximately by the corresponding 2x2 
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with dispersion relation 
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where I — ±1. The energy spectrum of this Hamiltonian 
is shown in Fig. Hfc) by the dashed curve making clear 
that this approximation is only valid for very small en- 
ergy (i.e. E/ji < 0.2) and near the Dirac point. The use 
of the approximate two band Hamiltonian allows for the 
extension of concepts as pseudospin to the trilayer sys- 
tem, which was originally defined for monolayer graphene 
(MLG) 2 . This extension makes it possible to derive sev- 
eral electronic properties from the conservation of pseu- 
dospin analogous to monolayer graphene^. 
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II. THE MODEL 



A. ABC Trilayer 



For ABC trilayer, the effective Hamiltonian near the 
Dirac point in one valley can be calculated using the tight 
binding formalism. Considering only nearest neighbor 
interlayer transitions, marked in yellow in Fig. QJa), one 
obtains^ 
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with a = (a x , Oy ) a vector of Pauli matrices, vf the Fermi 
velocity in monolayer graphene^ (vp ~ 1-01 x 10 6 m/s), 
k the wave vector and r describes the interlayer coupling 



Following a similar approach as before, the effective 
Hamiltonian obtained by a tight-binding model consider- 
ing only nearest neighbor interaction of the ABA trilayer 



is&: 



Haba + A = Hvf 



(f-k + S'I 2 t 

t* a-k r* 
r a-k-5'h 



in the same basis of orbital eigenfunctions as defined in 
Eq. I 2 is the 2x2 unit matrix and the 6' — 5/Hvp 
term corresponds to an externally induced interlayer po- 
tential difference of 6 which is described by the term A 
at the left hand side of the equation. This is a 6 x 6 
diagonal matrix given by 



A = Diag [5, 6,0,0, -S, -6] 



(7) 



Despite the strong resemblance with the ABC Hamilto- 
nian the ABA system is mirror symmetric with respect 
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to the central layer. Therefore, a unitary transformation 
that combines the orbital eigenfunctions symmetrically 
and antisymmetrically transforms the Hamiltonian into 
a block diagonal form: 25 



H' ABA + A = hv F 



a ■ k S'I 2 
5'I 2 o ■ k V2t 
V2t^ a ■ k 



(8) 



This new form of the Hamiltonian consists of a 2 x 2 
monolayer-like (top, left) and a 4 x 4 bilayer-like (bot- 
tom, right) block that are connected by the parts respon- 
sible for the interlayer potential difference. When <5 = 
the two blocks result in a superimposed linear (from the 
monolayer) and a hyperbolic (from the bilayer part) spec- 
trum near the Dirac point as shown by the solid curves 
in Fig. [TJd). In that case, electrons propagating in ABA 
TLG can propagate through two different modes, one 
monolayer-like and one bilayer-like mode. As long as the 
mirror symmetry remains intact, both modes will not in- 
teract and scattering between them is prohibited. When 
mirror symmetry is broken, e.g. by applying a different 
potential to every layer described by the term A, inter- 
band scattering is possible. In Fig. [TJd) the spectrum 
with non zero 5 is shown by dashed curves. Due to sym- 
metry breaking, the band crossing near the Dirac point 
is lifted making the linear monolayer-like bands become 
hyperbolic and the bilayer-like bands form a Dirac cone 
at the Dirac point while again crossing at higher wave 
vector. 



III. ELECTROSTATIC POTENTIAL 

In this paper we have considered two kinds of gates 
that influence the local electrostatic potential experi- 
enced by the electrons. The first one consists of a single 
gated device (SG) that causes a potential shift Vq equal 
for all three layers. The second one is a double gated de- 
vice (DG) that influences each layer differently inducing 
an interlayer potential difference 6 between neighboring 
layers. Both systems are translational invariant in the y 
direction. 

The SG and DG act as a boundary for which we calcu- 
late the transmission and reflection probabilities defined 
in the previous section as function of the angle of inci- 
dence and the Fermi energy of the incident electron for 
different configurations of the devices. If the transmitted 
electrons are measured inside the gated region, it can be 
modeled as a single boundary corresponding to a pn junc- 
tion. The ID potential profile is for convenience modeled 
by a step function 



Vpn 



(X) 



if x < Region I 
V I 6 + A if < x Region II 



(9) 



where Vq is the height of the potential and corresponds to 
the SG term and the term A in this potential describes 
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FIG. 2: (Colour online) Illustration of some of the different 
transmission and reflection channels for a generic two band 
system with an i th and a j th band. The potential and inter- 
layer bias in region II are Vo and S respectively, while both 
are zero in region I. 



the effect of a DG by inducing a potential difference be- 
tween the layers. This term is defined in Eq. (JT)), so the 
potential is Vo for the middle layer and Vo ± 5 for the 
top and bottom layer. The pn junction described here 
corresponds to the depicted schematic profile in Fig. [3] 
The single valley approximation used in this paper as- 
sumes that the potential varies over a length scale larger 
than the in-plane interatomic distance- a = 0.l42nm, 
but smaller than the electron wavelength. 

When the electrons are measured outside the gated 
region, one can describe it as a pnp junction or a potential 
barrier by 







^pnp (#) 



V I 6 



if x < Region I 
if < x < d Region II 
if x > d Region III 



(10) 



where d is the width of the gated region. 



IV. EIGENSTATES, CURRENT DENSITY AND 
TRANSMISSION PROBABILITY 

The eigenstates of these 6x6 Hamiltonians are six- 
component spinors consisting of a superposition of three 
times two oppositely propagating or evanescent waves 
characterized by three distinct wave vectors which we 
call fci , &2 and k$ . For a system that is translational in- 
variant in the y direction, the energy and k y dependence 
of these wave vectors h can be found from 



det [H(k x ,k y )-EI 6 }=0, 



(11) 



which leads to a sixth power polynomial in k x . The solu- 
tion of it corresponds to the inversion of the energy spec- 
trum. Since they solve the Dirac equation if 4" = E^f, 
the eigenstate spinors can be written as a product of ma- 
trices 



$>{x,y)=V£ (x,y)C, 



(12) 
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where V is a 6 x 6 matrix expressing the relative impor- 
tance of the different components of the spinor that can 
be constructed by solving the Dirac equation and 

c n- T ik-tx —ik-ix ikox — ikox ik-^x — ik-*x~\ —ik~.y 

t = Diag le , e , e , e 2 ,e 1 ,e 1 \ e yW , 

(13) 

where due to the translational symmetry in the y direc- 
tion, the y dependency is incorporated in an exponential 
phase factor and will be ignored from this point on. We 
denote the six component vector C as 



— ' ^ ^2 ' ^2 " ^3 " ^3 



(14) 



where the subscript i refers to the corresponding wave 
vector and the superscript +/— indicates the right/left 
propagating or evanescent states. The boundary condi- 
tions of the system under consideration will determine 
which of the components of the vector C are zero. 

Using the continuity equation, one can derive the cur- 
rent density for a general graphene multilayer Hamilto- 
nian with only nearest neighbor interlayer interactions 
and express it as a product of the above defined ma- 
trices. Such an Hamiltonian can be written in position 
representation as 



Hn 



(15) 



where a x ^ is a block diagonal matrix with n Pauli ma- 
trices <J x (y) on the diagonal and V consists of the other 
elements of the Hamiltonian that interconnect different 
atomic orbitals or induce a layer specific potential. Since 
the probability density is given by p = one can 

use the time dependent Dirac equation ihdt^ — H^> to 
obtain the current density: 

ihdtp = ih * + * f (<9t#)] (16) 

-v F (vtft) .a-r+*f) ^ 



(17) 



(18) 



where use has been made of the hcrmiticity of T and the 
a matrices. The current density is hereby given by 



] = V F 



(19) 



Introducing the matrix notation of Eq. (|12p . this expres- 
sion becomes 



where Ajj denotes the upper left element of the j th 
block. 

In this paper we consider the transmission of electrons 
between regions of different electrostatic potential or in- 
terlayer potential difference. This introduces a spatially 
varying potential term in the Hamiltonian in Eq. pip 
that determines those regions. The boundary is paral- 
lel with the y axis and we impose conservation of the 
transverse wave vector k y . This reduces the system to a 
one dimensional problem with the conservation of Fermi 
energy E and transverse momentum p y — hk y , which 
defines the angle of incidence on the boundary. We are 
interested in the dependency of the transmission and re- 
flection on these properties. 

Suppose the electron is incident at the left side of 
the boundary (region I) propagating in the k\ mode. 
The electron will be (partly) transmitted to the different 
modes at the left side of the boundary (region II) and 
(partly) reflected to the left propagating modes in region 

I. The boundary conditions at ±oo then yield expressions 
for the vectors C/ and Cu of the system in region I and 

II, respectively, as 



Cj = [l,rl,0,rl,0,rl] : 



(22) 
(23) 



In these expressions the coefficients V indicate the left 
propagating reflected waves and T indicates the right 
propagating transmitted waves. The subscripts denote 
the channel in which the waves are propagating and the 
superscripts indicate the incident mode of the wave. In 
Fig. [2] different possible channels are shown schemat- 
ically for a generic two band system. For TLG a third 
band needs to be included, leading to extra combinations 
of scattered transmission and reflection amplitudes. Us- 
ing the expression for the current density in the x direc- 
tion given by Eq. (|21|) , conservation of probability cur- 
rent leads to the normalization condition for the different 
channels as 



A x — Ir 1 
■^•1,1 |'i 

ki| 2 ^,i 



' A x 
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(24) 
(25) 



which allows to define the scattered transmission and re- 
flection probabilities as 



t'; 



AX 

A x 



and m 



A' 



(26) 



J 



J=v F C^£^V^aV£C. 



(20) 



Due to the properties of the spinor wave function, the 
matrix A = V^aP is diagonal consisting of traceless 2x2 
blocks that each correspond to a propagation mode. The 
resulting current density is 



A 



■3,3' 



(21) 



The transmission and reflection probabilities therefore 
depend on the value of the coefficients of the vector C. 
These coefficients can be found by matching the plane 
wave solutions of both regions at the boundary. Using 
the transfermatrix approach as explained by Barbier el 
al*^, one can model a sequence of different regions to 
create more complex structures such as a pnp junction. 

To find empirically relevant quantities, the transmis- 
sion probabilities can be used to calculate the con- 
ductance which is defined by the Landauer-Buttiker 
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formula^: 

j poo 3 

G(E) = G ^ / dk y J2 TL(E,k v ), (27) 

Z7T J -°° l,m=l 

with Go = 4e 2 //i, four times the quantum of conductance 
due to valley and spin degeneracy and L y is the length 
of the sample in the y— direction. 

Note that if both Dirac points are equivalent, the sys- 
tem is time reversal invariant. Therefore electrons that 
scatter reflectively from the i th to the j th band are equiv- 
alent to electrons reflecting from the j th in the i th band 
near the other Dirac point. A similar symmetry arises for 
the transmission when the potential, i.e. the eigenfunc- 
tions, of the first and last region of a series of boundaries 
are the same, such as with a pnp junction. Transmission 
from the i th into the j th band near the first Dirac point 
is then equivalent to transmission scattering from the j th 
into the i th band near the other Dirac point while being 
incident on the opposite side of the pnp junction. These 
equivalences, together with the equivalence between the 
Dirac points lead to the following symmetry in the reflec- 
tion and transmission probability 

R) = R{ and T] = T/', (28) 

where the second equation only holds if the first and the 
last region are equivalent. This decreases the number of 
different probabilities to only 6 reflection and 9 transmis- 
sion probabilities for a general potential and to 6 trans- 
mission and 6 reflection probabilities for a pnp junction. 

A final remark considering the symmetry of the system 
with respect to the sign of the angle of incidence has to 
be made. The Hamiltonians given in Eqs. ([I} and © 
are not symmetric under the change of the sign of k y . Al- 
though the electronic spectrum of these systems do have 
this symmetry, the obtained transmission and reflection 
results do not have to bare the reflection symmetry with 
respect to normal incidence. ABC TLG however has an 
additional symmetry, namely that by turning the system 
upside down and rotating in plane by an angle of 7r, the 
system is the same, but the sublattices ot\ and /?3, /?i and 
«3 and «2 and /?2 are interchanged. This transformation 
should not change the result, but it transforms the ABC 
Hamiltonian in such a way that k y — > — k y implying that 
the obtained results should also have this symmetry. A 
similar argument also holds for MLG, where both sublat- 
tices are interchanged, and for BLG, where the exchange 
is between the sublattices connected by the interlayer 
hopping and those that are not connected. Therefore, 
the symmetry also holds for unbiased ABA TLG since 
its Hamiltonian consists of MLG and BLG like blocks as 
shown in Eq. ([5]). These symmetry arguments are not 
valid anymore if a DG is applied since it breaks the in- 
terlayer sublattice equivalence of ABC TLG and the sep- 
arable behaviour of ABA TLG electrons into MLG and 
BLG like ones. Non symmetric results are thus expected 
for such systems. Although being counterintuitive, the 
occurrence of this asymmetry is not an unphysical result 
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FIG. 3: (Colour online) Schematic representation of the en- 
ergy spectrum of (a) ABC TLG and (b) ABA TLG at both 
sides of a SG potential boundary of height Vo without inter- 
layer bias. The colour of the curves indicate the wave vector 
associated with it. The minimum of the upper two bands of 
ABC TLG is located at energy a ~ 346meV . 

since the effect is exactly opposite for electrons near the 
other Dirac point. In the second valley, the electrons are 
described by a Hamiltonian that is similar to that near 
the first Dirac point, but with the exchange of E — > —E 
and k y — > —k y . The results for electrons (holes) in the 
first valley are therefore the same as for holes (electrons) 
with opposite transverse momentum in the second val- 
ley. The electron-hole symmetry of the system finally 
completes the argument. 

V. NUMERICAL RESULTS 

In the following sections we present results obtained 
for both SG and DG applied to trilayer graphene sam- 
ples with both ABC and ABA stacking configurations for 
both pn and pnp junctions. 

A. Single gated device 

ABC TLG has a maximum of three distinct modes of 
propagation, two of which are only propagating if the 

Fermi energy is large enough, i.e. E > a — j^J^Ji ~ 

346meV. In Fig. [H^a) the energy spectrum of ABC TLG 
is shown at both sides of the boundary. The branches of 
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FIG. 4: (Colour online) Transmission probability as function 
of the energy and transverse wave vector through (left) a SG 
pn junction of height V = O.571 and (right) a SG pnp junction 
of width d = 25nm and the same height using the two band 
Hamiltonian. 

the spectrum are coloured corresponding to the modes of 
propagation. 

At low energy, there is only one mode of propagation 
in ABC TLG which allowed before to introduce the two 
band Hamiltonian, Eq. that approximates the sys- 
tem as one with a cubic energy momentum relation. In 
Fig. [4] we show the transmission probability through a 
SG pn and a pnp junction for low energy in case of the 
two band system. Notice that at normal incidence, i.e. 
k y = 0, the transmission equals unity independent of 
energy or width of the barrier. This Klein tunneling is 
the consequence of conservation of pseudospin and oc- 
curs for rhombohedrally stacked multilayers with an odd 
number of layers. Furthermore, there is a region of k y 
corresponding to an angle of incidence of ±0 = 7r/6, for 
which the transmission is suppressed when E < Vb . This 
is another consequence of the pseudospinorial nature of 
the electrons in ABC TLG. At this angle of incidence, the 
propagating states inside the junction are disconnected 
from those outside resulting in a lower transmission. 

At non normal incidence for E < V, very narrow reso- 
nances show up for the pnp junction (see Fig. fj} that are 
similar to the Fabry-Perot resonances observed in MLG 28 
and BLG 1 ^-. The number of resonances and their energy 
and k y dependence vary with the width of the barrier. 
For the pn junction, spots of high transmission are found 
at non normal incidence. These spots are also a con- 
sequence of the chiral nature of the charge carriers and 
occur near an incident angle of cj) = ir/3£^ 

When the Fermi energy E of the electrons under con- 
sideration is larger than a, i.e. the minimum of the sec- 
ond band in Fig. the second and third modes of 
propagation become accessible. For these energies, the 
two band approximation is not sufficient anymore and 
interband scattering between the three bands needs to 
be taken into account. 

In Fig. [5]we show the transmission and reflection chan- 
nels for a pn junction of height Vq = I.571 as function of 
the transverse wave vector and the Fermi energy of the 
incident electron using the six band Hamiltonian. Note 



that Klein tunneling and cloaking when described by the 
two band approximation reoccur in the energy interval 
where only one band is propagating in both regions. This 
can be seen in the T-j 1 and R\ channel in the energy in- 
terval Vb — a < E < a. Outside this range, for lower 
energy the ki and k^ states inside the junction are prop- 
agating, giving rise to non zero scattered transmission 
probabilities T\ and T\ which reduces the transmission 
via the k\ channel. For larger energy, it is possible to 
reflect in the left propagating /c2 and k% states in region I 
which further reduces the direct transmission. For Fermi 
energy larger than the junction's height, the channel 
coincides with the two band calculation. Notice however 
that electrons impinging in the second or third band, are 
reflectively scattered in the third and second band re- 
spectively, rather than propagating in the k\ channel. 

In Fig. [6] the twelve transmission and reflection prob- 
abilities for electrons incident on a SG pnp junction of 
width d = 2bnm and height V = I.571 using the six 
band Hamiltonian are shown. The results for the 
and R\ channels are similar to those obtained from the 
two band Hamiltonian. Notice again Klein tunneling at 
normal incidence and the suppression due to cloaking at 
non normal incidence. Outside the energy range of va- 
lidity for the two band approximation, i.e. E < V — a 
and E > a, the structure of the resonances at low en- 
ergy is a superposition of the one similar to the two band 
system and another type of resonances. The latter are 
resonances due to propagation via the ki and k^ bands 
inside the junction region and therefore only show up 
when E < Vo — a. Furthermore, reflective scattering is 
large when all three modes are evanescent in the junction 
region. This is analogous to the interband scattering for 
BLGJ£ Another striking fact is that at normal incidence 
when a second mode of propagation is possible, Klein 
tunneling is suppressed in favor of scattered reflection in 
the R\ and R\ channels. This leads to the argument 
that Klein tunneling can be seen as the consequence of a 
suppression of the ability to backscatter rather than the 
link between forward propagating inside and outside the 
potential barrier. 

The large scattered reflection between the k^ and k% 
channels shown in Fig. [6] at normal incidence is a con- 
sequence of the tight relation between the ki and k% 
propagation modes. When Eq. (fTTj) is solved for nor- 
mal incidence, one finds an expression consisting of three 
hyperbolic bands, two of which intersect at E = 71. At 
normal incidence, the reflected branch of the £3 spectrum 
as indicated in Fig. 02a) therefore corresponds with the 
forward propagating ki band and vice versa. Although 
at non normal incidence this relation is no longer valid, 
the reflective scattering remains strong at near normal 
incidence. 

In Fig. [7f a) we show the energy dependence of the con- 
ductance for the same system as for the previous results. 
The resonances that are visible in the transmission prob- 
ability show up as peaks in the conductance. Further- 
more, the results indicate the availability of additional 
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FIG. 5: (Colour online) Transmission and reflection probabilities for a single gated pn junction on ABC TLG of height Vo = 1.671 
as function of the energy and transverse momentum k y . 



propagation modes via the and channels for higher 
energy resulting in a clear signature of higher conduc- 
tance. The increased conductance for the high barrier is 
a distinct feature from the results of the two band ap- 
proximation. For low energy (E < O.571) and just below 
the barrier's height (71 < E < V), the conductance is 
raised due to propagation via the second and third chan- 
nels in the barrier region. This is absent in the two band 
approximation. 

The SG imposes a potential that treats each layer in 
the same way and thus keeps the existing symmetries 
of the system. As a consequence, the spectrum of the 
system is only shifted by an overall potential term as 
shown in Fig. EJb). The monolayer-like and bilayer-like 
propagation in ABA TLG therefore remains the same 
and so the electrons are described as if they propagate in 
MLG or BLG. The transmission and reflection through 
pn and pnp SG devices on MLG and BLG has already 
been investigated in depth 2 - ' 16 i 28 . It was pointed out that 
the MLG electrons exhibit Klein tunneling when they 
hit the boundary perpendicularly. They are transmitted 
with unit probability, irrespective of the height \ or the 
width d of the gated region. However, BLG electrons are 
cloaked from the propagating states at this angle of in- 
cidence and therefore their transmission is suppressed 15 . 
At non normal incidence, cloaking and Klein tunneling 
can occur if certain conditions are satisfies as described 
recently 2 ^. The MLG and BLG like propagation channels 



have each their own separate contribution to the conduc- 
tance of an ABA TLG sample. In Fig. [TJb) we show the 
energy dependency of the conductance of an ABA TLG 
sample with a SG pn and pnp junction. From the con- 
tributions of the separate channels one finds that when 
E < Vo, the MLG conductance G\ is larger than that of 
the BLG due to the Klein tunneling, but the BLG con- 
ductance G\, although being cloaked, shows clear peaks 
reminiscent from the BLG transmission resonances. For 
larger energy, when E > Vo, the BLG contribution is 
however larger than that of the MLG. The overall con- 
ductance of a pn junction doesn't differ much from the 
pnp junction when E > Vo, while it misses the resonances 
at energy lower than the junction's height. 



B. Double gated device 

A potential difference between the layers of an ABC 
TLG destroys the degeneracy of the valence and con- 
duction band at the Dirac point by opening up a band 
gap between the m 10 ' 29 . This is illustrated in Fig. [5Ja). 
Due to the absence of propagating states inside the band 
gap, the transmission is reduced in all channels which 
can be seen by the low conductance as shown in Fig. 
[H£a) where both the pn DG and pnp DG conductance is 
shown. When the propagating states become available, 
the conductance increases sharply to unity in the regime 
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FIG. 6: (Colour online) Transmission and reflection probabilities for a single gated pnp junction on ABC TLG of height 
Vo — I.571 and width d = 25nm as function of the energy and transverse momentum k y . 



where only one band is available to conduct. As more 
bands become possible to conduct, additional transmis- 
sion channels contribute to the overall conductance in- 
creasing it to almost perfect conductance. Due to the 
absence of resonances, the pn DG and pnp DG conduc- 
tance are similar, only showing the effect of additional 
transmission channels. The DG lifts the valley degener- 
acy and therefore the symmetry in the scattered trans- 
mission expressed in Eq. (I28[) is no longer valid but are 
the same up to a reflection with respect to normal inci- 
dence. The conductance however sums over both posi- 
tive and negative angles so the effect is absent. Since the 
asymmetric effect is only possible when several modes 
of propagation are available, it will be more pronounced 
when considering ABA TLG. 

The breaking of the interlayer symmetry induced by 
a DG couples the MLG and BLG modes as shown in 
Eq. ([8]). We have calculated the scattering between the 
linear and the lowest hyperbolic band due to a pn DG 
and a pnp DG. The way in which the ABA spectrum 
is influenced by the DG is shown schematically in Fig. 
IHJb). A band gap is created between the linear bands 
for E < 5, the hyperbolic bands however remain gapless 



touching linearly at the Dirac point and intersecting in 
a circle of radius k = 5/jia around the Dirac point. For 
low energy, this feature allows for three conduction modes 
of the same band while for a Fermi energy that is a little 
higher, only one mode is available. 

In Fig. [TOT a) we show the angular and energy depen- 
dency of the transmission and reflection probabilities at 
a pnp DG. They can be divided in several regions defined 
by which mode is propagating as shown by the dashed 
lines superimposed on the results. A clear feature is the 
hyperbola reminiscent of the gapped Dirac cone due to 
the interlayer asymmetry in the junction region which 
is displayed as a white dashed curve. Inside this hyper- 
bola, it is possible to propagate through the ki channel 
and this increases the linear to linear energy band T\ 
transmission. In contrast to MLG electrons, it is pos- 
sible to be backscattered in the linear band. A second 
clear feature is the Dirac cone defining the MLG propa- 
gating states outside the junction, which is indicated by 
a black dashed curve. Outside the Dirac cone, only the 
BLG-like states are propagating both inside and outside 
the junction. This leads to angular depending resonances 
that vary with the length of the pn junction. The previ- 
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FIG. 7: (Colour online) Energy dependency of the conduc- 
tance at a SG junction for (a) ABC TLG with potential 
height Vo = 1.671 and (b) ABA TLG with potential height 
Vo = O.571. The solid black curves correspond to the to- 
tal conductance through a pnp junction of length d = 25nm, 
the dashed coloured curves indicate the contributions of the 
different propagation channels. The dotted black curve corre- 
sponds to the conductance through a pn junction of the same 
height. 



ously discussed angular asymmetry is clearly seen in the 
reflection channels and the scattered transmission chan- 
nels. While the scattered reflection keeps the interband 
symmetry R\ = R\, the scattered transmission differs. 
In Fig. rTUT b) we show the difference between both scat- 
tered reflection channels. The result makes clear that the 
scattered transmission is the same under a flip of the sign 
of k y , i.e. T^iky) = Tl{~k y ). This is an immediate con- 
sequence of the time reversal symmetry of the system. 
As mentioned above, an electron scattering from band 1 
to band 2, near the Dirac point K, T 2 X K , is equivalent 
to one scattering from band 2 to band 1 near the other 
Dirac point K', T^ K ,. Since the latter is the same as 
Ijjj- with a sign flip in k y , the scattered transmissions 
near the same Dirac point are the same under a change 
of sign in k y . Note that these asymmetric results are not 
carried through in the conductance. 

In Fig. IHIb) we show the conductance for a pn and a fi- 
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FIG. 8: (Colour online) Schematic representation of the en- 
ergy spectrum of (a) ABC TLG and (b) ABA TLG at both 
sides of a DC boundary with interlayer bias S. The colour of 
the curves indicate the wave vector associated with it. 



nite width pnp DG in ABA TLG. The resonances caused 
by the T| channels show up as small bumps in the con- 
ductance when the Fermi energy is below S. For higher 
energy, the contribution of the MLG like mode of propa- 
gation is marked by a rise in conductance due to the con- 
tribution of the G\. In Fig. [IT] the conductance is shown 
as function of the width of the pnp DG with strength 
5 = O.571 for different levels of the Fermi energy. In 
Fig. ITTT a). we show the conductance for a Fermi energy 
E < S. Now there is only one mode of propagation in the 
junction and the oscillatory behaviour is only due to the 
resonances of the G\ term. The direct G\ contribution 
diminishes exponentially since and obtains a finite value 
independent of the width of the junction, similar to the 
interband scattering contributions. Fig. [TTTb) shows the 
conductance for E > 8. Now the contributions of both 
direct channels, G\ and G§> are almost equal and oscil- 
late with the width of the junction while the scattered 
transmission oscillates oppositely. In this way, one can 
determine the amount of electrons scattered between the 
monolayer- and bilayer-like bands by applying a DG of 
the correct strength and width. 
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VI. CONCLUSION 

In this paper we have studied the electronic transport 
in trilayer graphene through pn and pnp junctions con- 
sisting of a single and a double gate. For ABC TLG 
we have shown that the availability of different modes 
of propagation at higher energy leads to typical features 
in the conductance and that Klein tunneling occurs only 
if there is no other propagation mode available. Fur- 
thermore, we have calculated the effect of the band gap 
induced by a DG on ABC TLG in the conductance and 
shown that the conductance is nearly zero in the gapped 
region. 

We also modeled the effect of SG and DG gates on 
ABA TLG and found that the SG behaviour is a super- 
position of a monolayer and bilayer like system. The DG 
however mixes both types of bands and breaks the an- 
gular symmetry with respect to normal incidence. This 
peculiar result emphasis the necessity to include both 
Dirac points even if intervalley scattering is prohibited 
since the electron behaviour near the other Dirac point 
restores the symmetry. The DG finally allows electrons 
to scatter between the monolayer and bilayer like bands 
and we have shown that the strength of scattering de- 
pends on the with the width of the pnp junction and the 
Fermi energy considered. 



FIG. 9: (Colour online) Energy dependence of the conduc- 
tance at a DG junction of strength S = O.371 for (a) ABC 
TLG and (b) ABA TLG. The solid black curves correspond 
to a pnp junction of length d — 25nm, the dashed coloured 
curves indicate the contributions of the different propagation 
channels. The dotted black curve corresponds to the conduc- 
tance through a pn junction of the same strength. 
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FIG. 11: (Colour online) Conductance as function of the 
width of the pnp DG of strength 8 — O.571 for Fermi energy 
(a) E = O.371 and (b) E = O.871. 



